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o ■ 

^ ' Abstract 

(N : 

^ , Vector-Meson-Dominance model contribution to vr^ 47 is calculated. The 

result confirms old estimates that this contribution is much smaller than the purely 
' electromagnetic photon spliting graph contribution calculated earlier. 

a^ : PACS: 12.40.Vv, 13.25. Cq 

^ ! In the framework of the Vector-Meson-Dominance (VMD) model, the decay vr^ 47 

Q\ [ was considered in [Tj and 0. However, no complete evaluation of the decay width can be 

^ I found in these papers, because the calculations are too tedious to be performed by hand. 

^ ■ Only an upper limit R < (7-i-8.6) x 10~^^ was given in |T| as a result of partial evaluation 

^ ■ of the squared decay amplitude. Here R stands for the ratio 

o; _ r(7r° ^ 47) 

r(7r0^27)' 

OhI In |2] the completely different number R ~ 10~^ is quoted, but we believe this result is 

^ ' erroneous. 



p. 
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The present experimental upper limit 2 x 10^^ on the tt^ 47 decay branching 
ratio was obtained long ago 0. Recently a Letter of Intent for a new experiment at 
PSI appeared [5 aimed to search C-noninvariant decay 7r° 87. A by-product of this 
5^ I experiment will be an improved measurement of the allowed decay 7r° 47 as the most 

important background jl]. Taking in mind this experimental situation, the complete 
calculation of the VMD contribution seems desirable. 

Using the advantage of the REDUCE Computer Algebra System 0, we were able to 
perform this complete calculation and the results will be presented below. 
For the standard tt^ 27 amplitude (/t^ ~ 93 MeV) 



the corresponding decay width looks like 

r(7r° ^ 27) 



2 S 



327r3/2 



While, by using the Kumar's parametrization of the covariant phase-space jH] , the 7r° 47 
decay width can be written in the form jHl E] 

r(7r° ^ 47) = 



1 



4! 2iV«mJ i J J JA(m2,32,s'2) -f- •', Vl - C" 

1/ V 



where various Mandelstam-type variables are defined in Kumar's paper [Oj. 
In the VMD model, the 7r° 47 decay amplitude has the form jl] 

M = r(l, 2; 3, 4) + T(2, 1; 3, 4) + T(3, 2; 1, 4) + T{2, 3; 1, 4) + 

+r(4, 2; 3, 1) + T(2, 4; 3, 1) + T(l, 3; 2, 4) + r(3, 1; 2, 4) + 

+T(1, 4; 3, 2) + T(4, 1; 3, 2) + T(3, 4; 1, 2) + T(4, 3; 1, 2), 

where for our choice of the coupling constants (note that V ^ tt'J couphng constant is 
defined as egv-K-r) and up to an irrelevant phase 



T(l,2;3,4) 



a 



Si — rn?p Si — 



S2 - ml 



X 



^^^lv\a4K^AK (^ais-ysih + h + k^f elk{ eafs'j'S'ih + ^3 + kiY'etk2 . 

While calculating the total decay rate, the permutation symmetry of the phase space can 
be used to classify the 144 terms in the squared amplitude into 7 types [T] : 

• |T(1,2;3,4)|2 - 12 diagonal terms, (i) 

• r(l, 2; 3, 4)T(2, 1; 3, 4) - 12 terms, (n) 

• T(l, 2; 3, 4)T(1, 3; 2, 4) - 24 terms, (iii) 

• T(l, 2; 3, 4)T(2, 3; 1, 4) - 24 terms, (iv) 

• T(l, 2; 3, 4)T(3, 1; 2, 4) - 24 terms, (v) 

• T(l, 2; 3, 4)T(3, 2; 1, 4) - 24 terms, (vi) 

• T(l, 2; 3, 4)T(3, 4; 1, 2) - 24 terms. (vn) 

All terms of a given symmetry type contribute equally in the total decay rate. Terms of 
the type (vii) turn out to be zero after doing the photons polarization sums. 

Introducing a dimensionless version of Kumar's invariant variables and using the same 
notations for them: 



1 1 

si = -^{q - kif, S2 = -^{q - ki ~ k 
mi m 



111 

ui = -^{q - /C2)^ U2 = -^{q - ksf, t2 = -^{q - k2 - k^f, 
m^ m^ m^ 

we get after performing the polarization sums by REDUCE 

r(7r° ^ 47) a 



r(7rO 27) 4(47r)' 



■X 
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1 Si 1—S1+S2 



X 



/ dsi / dS2 / , ^ = 



F{si,S2,Ui,U2,t2{C))- 



Note that the invariant variable t2 is a hnear function of the integration variable 
where 



rj 



X{l,S2,s',)-{l-Siy + {l-U,y 

2(l-iii)^A(l,S2,4) 

(l-4)2-(l-«2)^-A(l,.2,4) 



2(1-M2)VA(1,S2,4) 

and A(a;, y, z) = x"^ + + z"^ — 2{xy + xz + yz) is a conventional triangle function. Besides 

s'2 — 1 -\- S2 — Ui — Si, s'^ — 2 — Si — Ui — U2- 
The hmits of integration for the ^2- variable are 

uf ^1- ^{ui + Si) ± ^^A(1,S2,4)- 

The function S2, mi, M2, ^2) can be decomposed into six parts, each of them corre- 

sponding to the particular symmetry type mentioned above: 



F = Fi + F2 + F3 + - [F4 + F5 + Fe] . 



Here 



+ 



a,„ 1 



si-bp si-bu,\ (s2 - 1)^ 



Pi, 



F4^ 



si - bp si- b^ 



+ 



dp ^ dfjj 



ui - bp ui- b^ 



{S2 - If 



Si - bp 61 - k 



Pz, 

S2-I Si- S2+Ui+U2~t2-2 



Si - bp Si - 



F.^ 



+ 



ui — bp ui — k 



1 



si - bp si- b^ 
1 



+ 



S2-I t2-l 



Pa. 



U2 - bp U2- b 



X 



X 



S2 - 1 Si - S2 + Ui + U2 - t2 - 2 



dp ^ dijj 



Si — bp si — b^ 



+ 



U2 — bp U2 — b 



S2-I t2-l 



Pe, 



where 



bp = my ml, b^ = ml /ml, 



and Pi, i = 1 -h 6 are certain polynomials of the variables Si, S2, Mi, ^2 given in the 
appendix. 

For numerical calculations we need gv^K-y coupling constants and they can be estimated 
from the V tt^j decay widths with the following result [H]: 



ap ^ 0.0097, ^ 0.098, bp ^ 32.53, ^ 33.55. 
After the numerical calculations, we obtain 



in agreement with the estimates given in Therefore, the VMD model contribution in 
the 7r° — > 47 decay width is indeed very small, many orders of magnitude smaller than 
the photon splitting graph contribution R ^ 2.6 x 10~^^ [H]. 

Recently hadronic contribution to 7r° 47 was estimated by using constituent quark 
loop model ^0]. Earlier this contribution was calculated by using chiral perturbation 
theory ^1]. In both approaches the hadronic contribution was also found to be negligible. 
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Appendix 

The expressions for the Pj(si, S2, mi, ^2, ^2), i = I -^6, polynomials are: 

Pi = I sllsj + (1 - + 2(S2 - - Si - Ml)] + S2(S2 - 2MiSi)| 

-P2 = 52 I S2 + SiMi(si - 2s2 + 3mi - 2) + SiMi[(l - UiY + sj- 252^1]! 
^3 = stisj + l) + sl\ sli2 - 3ui -U2 + 2t2) - 2sl + Si(m2 - Ml - 2) + 



+2(1 - Ml - M2 + ^2) j + S2 I + sliAul +U2- 3t2 - 2) + 

+ sl{4ul +ul + tl- AUit2 + 3MiM2 - ^2^2 - 2mi - 2M2 -t2 + 2) + 
+ Si{2u\ -u\- 2Uit2 + U1U2 + ^2^2 + 3^2 - 2t2 - 2) + 

+ (1 - Ml - M2 + t2f I + S2S1 I s\{t2 - Ui) + 

+ s\{bUit2 - -t\- 2UiU2 - t2U2 + Ml + 2t2) + 

+ Si(l - Ml - M2 + t2)(3Mi - Mit2 + U1U2 + Mi - 3t2)- 



gp^^ ^ 0.73 GeV"\ g^^^ ^ 2.32 GeV 



which translates into 



P ^ 3.3 X 10 



16 
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-Ui{l - Ui- U2 + I + s\{u\ + SiMi - 51^2 - «1^2 + «1«2 - 



{si + + ^2)^ - 2^2(^2 + 2mi)- 



-251^2(1 -\-ul-\-u2- 12) 



I + 2s2lil I s\[t2 - U1U2 - {Ui - ^2)^] + 



+Sit2{u2 — 2t2 - 1) + SiUi 



Ui{3t2 — SU2 — 2Ui + 1) +U2 + 2t2 — {U2 — ^2)' 



-Ui{l - Ui- U2 + t2){2t2 - Ml - U1U2) I + 



+Ui{si + Ui) 



Slt2 + - Ml - M2 - Si + ^2) 

P5 = 2S2(S1M2 + 1) - I 2slu2 + Si{ul + 3UiU2 - 2u2t2 + Ml - 3x^2 + 2)- 

-2(1 - til - + ^2) I + S2 I 2s\u2+ 

+ sl{'u\ + + 8M1M2 - 6M2t2 + 4?7,i - 6?7,2 " 4t2 + 2) + 



+Sl 



1*2 + 2^2(3^1 - ^2 - 2) + U2{lul - 8Uit2 - 6ui + 2tl + At2 + 8)+ 



+A{ul - Uit2 -t2-l) 



+ 2(l-Mi-M2 + t2)^| + 



+2SiS2 \ sl[2t2 - Ml(l + Ml + M2 - ^2)] " Si 



-2Ui{l + U2) (1 - M2 + 2^2) + ^2(2 + M2i2 - 2M2 + 2^2) 
+ (1 —U1—U2 + ^2) 



+ 



+ 



Ui{2UiU2 — U2t2 + M2 + Ml - ^2 - 1) - M2i2 
+Si(si + M2)[Mi(1 - Iti - lt2 - Si + t2) + Sit2]^, 
Pe = S2 I Siti2[(Mi + ti2)^ - 2^2(1 + Ml + 1*2 - ^2)] + 2^2 - 2t2Ui{2Ui + ^2) + 

+ Mi(mi + M2)^| + 2S2 I M2S^(2Mit2 - M1M2 + ^2 - M^ - tl)- 



-Sl 



tl{uiU2 + Ml + U2) — t2{ul + 2u\u2 + 2u\ + 2x41^2 + '^2 ~ '^2) + 



+ Mi(wi + M2)(m^ + M2 + M1W2 - M2) 



+ 



+u\{l -U1-U2 + t2){u\ + U1U2 - 2^2) I + 
+ {u\ + SiW2)[si(i2 - Ml) + Ui{l — U1—U2 + t2)f'. 
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